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Abstract 

We consider the Heisenberg spin chain in the presence of integrable spin defects. Using 
the Bethe ansatz methodology, we extract the associated transmission amphtudes, that 
describe the interaction between the particle-hke excitations displayed by the models and 
the spin impurity. In the attractive regime of the XXZ model, we also derive the breather's 
transmission amplitude. We compare our findings with earlier relevant results in the context 
of the sine-Gordon model. 
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1 Introduction 

The presence of impurities in physical systems is an issue of great consequence, especially 
when dealing with more realistic models, and confronting experimental data. Along this 
spirit, integrability offers a framework where impurities may be naturally incorporated in a 
controllable manner. Although there have been numerous recent advances in both quantum 
[l]-[8], and classical [9]- [19] models with integrable defects, many questions still remain 
open. In the present investigation, we restrict our attention on quantum spin chains in 
the presence of a single integrable defect and extract the physical information concerning 
scattering processes within such models, directly from the Bethe ansatz equations. Our 
approach thus aims at complementing and enhancing the picture described in earlier works 
[4] in the context of integrable field theories. 

The algebraic frame describing the presence of a point-like defect in a discrete integrable 
theory is by now well established through the quantum inverse scattering method (QISM) 
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[20, 21]. The formulation is based on the existence of a defect Lax operator that satisfies the 
same quadratic quantum algebra as the bulk monodromy matrix. In general, let us consider 
a one dimensional (A^ + l)-site theory with a point like defect on the nth site. In this case 
the modified monodromy matrix of the theory reads as 

T(A) = Wi(A) Lo^(A)...Lo„(A-e)...Loi(A), (1.1) 

where L corresponds to the "bulk" theory, L corresponds to the defect and 6 is an arbitrary 
constant corresponding to the "rapidity" of the defect. Both Lax operators satisfy the same 
quadratic algebra 

i?i2(Ai - A2) Li(Ai) L2(A2) = L2(A2) Li(Ai) R^^iK - A2) , (1.2) 

where the i?-matrix is a solution of the Yang-Baxter equation (see e.g. [21] and references 
therein). The monodromy matrix of the theory T(A), naturally satisfies (1-2), guaranteeing 
the integrability of the model. The Hamiltonian of any generic system with a point-like 
defect is given by 

Af+l 

n^-[Y. ^n+M + L+in{0) L-l,^iO) + Ln+miO) k-m+m L^j^jO)) (1.3) 

the "dot" denotes the derivative with respect to the spectral parameter. We focus here on 
the situation where -^^(A) = -R(A), also define R = VR, V is the permutation operator. Recall 
that the R matrix reduces to the permutation operator at A = 0. 

Note that here we are going to focus on the ant i- ferromagnetic regime of the XXX and 
XXZ models. The derivation of the Bethe ansatz equations (BAE) is straightforward in 
the case where highest weight states exist. Amongst others, the thermodynamic limit of the 
BAE provides us with the scattering information for the given model, which we exploit in 
order to derive our results. Generalization of our results in the presence of multiple defects 
is straightforward within the QISM frame. 

In the subsequent sections, we investigate the interaction between the particle like exci- 
tations displayed in the XXX and XXZ spin chains and the defect. These interactions are 
described by generic transmission matrices, that satisfy the quadratic algebra (1.2) with an 
overall physical factor, which is explicitly computed by means of the BAE. The XXZ model 
with a defect is studied in both the attractive and the repulsive regime. In the latter the for- 
mation of bound states between solitons and anti-solitons, called breathers, is allowed. After 
describing the scattering process for the breathers, we derive the transmission amplitude 
between a breather and the defect of the theory. We also compare our findings with ear- 
lier results, reaching complete agreement. Finally, the appendices contain several technical 
points and physical checks that further confirm our results. 
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2 The isotropic case: XXX model 



We begin our analysis considering the isotropic XXX spin chain in the presence of a single 
defect (see also e.g. [22] and references therein). Let us first recall the quantum Lax operators 
for the bulk model, L, and the for the defect L. The generic defect matrix is given by 



L{\) 



A + + f 



iS- 



(2.1) 



where the algebraic objects S^, are the generators of the SU2 algebra, 



±s 



± 



s- 



2S' 



(2.2) 



The bulk L-matrix is actually the i?-matrix of the model, which corresponds to the spin-| 
representation of the SU2, that is the following identifications are implemented in (2.1): 



± 



(2.3) 



As usual, cr^, denote the familiar 2x2 Pauli matrices. In the finite case, for an n = 25* + 1 
dimensional representation of spin S, the algebraic objects and S*^ are represented by 
n X n matrices defined as 



k=l 



n-l 
k=l 



11-1 

Ck Ck+ik 

k=l 



(2.4) 



where we define the matrix elements: (en 



ijjkl 



5ik 5ji and 



^(n + l-2fc), Ck = ^/k{n - k) 



(2.5) 



This choice may be thought of as the isotropic analogue of the type-II defect studied in [4], 
which will be analyzed in the subsequent sections within the XXZ spin chain context. The 
findings of this section may be seen as the quantum discrete analogues of the results on 
the Landau-Lifshitz model [19]. Classical scattering in the context of the Landau-Lifshitz 
model should be considered and comparison with our findings should provide a more concrete 
correspondence between the classical and quantum models in the presence of defects. The 
Hamiltonian of the model is given by (1.3) where L(0) = I, 



Rjj+i{0) 



2 v^J^J+i + 44+1 + 44+1 + Ij%+i j , 



(2.6) 
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The generic defect matrix as well as the bulk L-matrix possess highest weight states 
{S'^ 1+) = 0), thus the typical algebraic Bethe variation may be applied and the correspond- 
ing Bethe ansatz equations (BAE) are immediately obtained^ 

N 

ey{\i - 6) ef (A,) = - J] e2(A. - A,), y = 2S, (2.7) 

i=i 

where is the rapidity associated to the defect, and we also define: 

A + — 

e.(A) = . (2.8) 

Having the Bethe ansatz equations at our disposal, we are now in a position to derive the 
physical quantities which describe the scattering processes on the chain. 



2.1 The transmission matrix 



Before we proceed with our analysis we shall recall that in the thermodynamic limit the 
solutions of the BAE may be expressed as "strings" with a real and an imaginary part. This 
is based on the so called "string hypothesis" , stating that the Bethe roots may be cast as 

A^'^'^) = Ao + ^(n + 1 - 2j), J e {!,..., n}. (2.9) 

Recall also that the total spin of a state may be obtained through the asymptotic behavior 
of the transfer matrix, and is given by the following familiar expression: 

S' = j + S-M. (2.10) 

The spin of the state should be non-negative, thus the restriction M < y + is manifest 
(below the "equator"), while the rest of the states can be obtained by starting from a 
reference state, which is a "lowest weight" state {S~ |— ) =0). The energy and momentum 
may be also explicitly expressed in terms of the Bethe roots {Aj} (for more details on the 
Bethe ansatz formulation and relevant physical expressions see also e.g. [20, 23, 24]). 

As is well known, the anti-ferromagnetic ground state is a "filled Fermi sea" with real 
solutions, i.e. 1-string configurations. It is easy to check that the ground state has an overall 
non zero spin = S, where S = S — ^ is the "shifted" spin, which naturally emerges 
through the Bethe ansatz approach, as will be transparent subsequently when deriving the 
transmission matrix. Implementation of suitable string configurations can "correct" the 

^The BAE are valid for any S ^ real number. 



4 



spin, i.e. provide a total spin zero, without modifying the state's energy. In general, it is 
straightforward to show that in the presence of an n-string the spin becomes 

S' = S-{n-l), (2.11) 

while the energy and the momentum of the configuration are left intact. 

A generic state contains m particle-like excitations, which are interpreted as m holes in 
the "filled Fermi sea" . The density of the state may be obtained then, here we follow the 
standard formulation [20, 25] for the derivation. Starting from the Bethe ansatz equations, 
and after taking the logarithm and differentiating, we define the density of the state in the 
presence of m holes in the filled Fermi 

(7(A) = ai(A) + -ay{X - 6) - / dX' a2(A - A') aiX') + a2(A - A,) , (2.12) 

where we have defined 

an(A) = ^^ln(e„(A)). (2.13) 

Passing to the thermodynamic limit, we have exploited the following basic formula in the 
presence of m holes, with associated rapidities A^: 

T7 E /(^^■) ^ / /(^) ^(^) - ^ E /(^^■) • (2-14) 

j=l J-oo 

The Fourier transform of is needed subsequently, and is given by^ 

a„(w)=e-"^. (2.16) 

For our purposes here, we shall focus on the case where m = 2 (A^ is assumed to be odd), 
in order to derive the "kink" scattering matrix, as well as the transmission matrix. Using 
the machinery described above, the density may be expressed in a compact form as 

a(A) = ao(A) + - ( E ^^(^ - ~^k) + n^X - 9)) , (2.17) 
fc=i 

with (To (A) being the density of the ground state. The Fourier transforms of the latter 
quantities have been computed 

1 -(y-i)^ 
^'^("^=2^^' '^^") = 2^;h(f)' = ^'-''^ 



^We have considered the fohowing conventions for the Fourier transformations: 

1 

2^ 



/(A) = ^ I dio e-^-^f{u:) , /(c.) = / d\ e'-V(A) • (2.15) 
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Recall also that 



ao(A) = e{X), 



and 



1 dp{X) 



(2.19) 



27r dX 

with £ and p being the energy and the momentum of the hole excitation (kink), respectively. 
Recall also that 



a{\) 



1 dh{X) 



(2.20) 



dX 

h{X) is the so-called counting function and h{Xi) = Ji, where Ji are integer numbers. 

In order to identify the scattering amplitude between two holes as well as the hole-defect 
transmission amplitude, we compare the expression providing the density of the state (2.17) 
with the so called quantization condition, with respect to the hole with rapidity Ai: 



5(Ai,A2,e)-l)|Ai, A2) = 



(2.21) 



where S = e**, and p{Xi) is the momentum of the respective hole. Comparison of the 
quantization condition with the state's density (2.17) would provide the "kink-kink" scat- 
tering amplitude as well as the transmission amplitude, given that the factorization of 
the scattering is evident (see also [26]). More precisely, the study of the one-hole state 
would simply provide the transmission amplitude, which physically describes the inter- 
action between the particle-like excitation and the defect, thus factorization of the type: 
S'(Ai, A2, 6) = S's(Ai,A2) T(Ai,9), in the case of the two-hole state is manifest. Keeping 
these considerations in mind, the kink-kink amplitude as well as the transmission amplitude 
for the XXX model with a single defect can be derived then as 



Ss{X) = exp 



du 



T(A) = exp 



du 



(2.22) 



where A = Ai — A2 and A = Ai — 9. 

Taking now into account the following useful identity 



1 

2./0 



duj e 2 



In 



r(^) 



u coshf '"r(^; 



(2.23) 



the familiar expression for the XXX kink-kink scattering amplitude [20, 27] is reproduced in 
terms of F-functions 

^ ' r(-f + 1) r(f + 1) ^ ' 

This is the first eigenvalue of the scattering matrix and it is three-fold degenerate due to the 
underlying 5U2 symmetry. There is one more eigenvalue corresponding to the singlet state, 
which may be derived by considering the state with two holes and one 2-string with a real 
center: 

Ai + A2 



An 



(2.25) 
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The second eigenvalue may be identified then, and turns out to be 

Ss{x) zA + r 

The S'-matrix, which satisfies the Yang-Baxter equation, may be then cast as: 

ftX + l \ 

iX 1 
1 iX 

V iX + lJ 



(2.26) 



(2.27) 



which is the famihar XXX scattering matrix [20]. Unitarity and crossing are also explicitly 
checked and are satisfied by the extracted S'-matrix (see Appendix A). 

The transmission factor may be also identified through expression (2.22), and is found 
to be given by 

T(A,^) = ^^- (2.28) 

r(f + f + i)r(-f + 1 + 1) 

The spin associated to this state is identified through the BAE, and is given as: 5*^ = 
5" + !, (5 = 5'— |), corresponding to the highest spin 5*^ eigenvalue (Appendix A). Notice 
the "shift" of the spin which naturally emerges through the Bethe ansatz process. We have 
been able so far to identify the first eigenvalue of the transmission matrix. As shown in the 
Appendix A, there are only two distinct (li = 25" + l)-fold degenerate eigenvalues associated 
to the generic transmission matrix. 

As already pointed out, in order to identify the transmission amplitude, it is sufficient to 
consider the state with one hole (A^ even) (see also e.g. [28]). The other distinct eigenvalue 
may be found using suitable string configurations that modify the spin, but leave the energy 
of the sate intact. Considering n-strings, suitably positioned with respect to Ai and 9, 
one may determine each eigenvalue corresponding to the appropriate spin eigenvalue 

G jiS^jn, ■ ■ ■ , S + 1/2|. The rest of the negative spins may be obtained starting from 
a reference state, which is a "lowest" weight state (see also Appendix A). For instance, it 
is easy to identify the spin eigenvalue of a state with one hole of rapidity Ai, and one 
ra-string: nG |2,...,no|, where we define 

~ 3 

no = S + - =^ S'^j„ = 0, S half-integer 

no = ^ + 1 ^ = 1, S integer. (2.29) 



The spin turns out then to be 



S' = S+^-{n-l). (2.30) 
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The real center of the n-string is positioned at 

Ao = ^Ai+®, x = ^. (2.31) 
X X n — I 

The position of the ra-string is determined using suitable quantum numbers, which char- 
acterize the state. We shall not give the technical details of the proof here, however the 
interested reader is referred to [24] for more details. Given the arguments above, the second 
eigenvalue of the transmission matrix can be derived for such a state: 

IMM>^'^-^-^ (2.32) 

T{\,S) i\ + S + l 

corresponding to the spin eigenvalue (2.30). 

Having the eigenvalues of the transmission matrix at hand, and bearing in mind that for 
purely transmitting defects the following quadratic algebra is satisfied [1] 

Si2(Ai - As) Ti(Ai) T2(A2) = T2(A2) Ti(Ai) §i2(Ai - A2) , (2.33) 

where S is given in (2.27), we conclude that 

This is the generic matrix associated this time to the spin S = S — ^ representation of SU2. 
The "shift" of the spin via the Bethe ansatz process, is once more pointed out. 

The findings of this paragraph can be further confirmed by checking the basic require- 
ments of unitarity and crossing symmetry, and of course by recalling that the transmission 
matrix has to satisfy (2.33). The transmission matrix presented above clearly satisfies the 
fundamental algebraic relation due to its structure. Unitarity and crossing symmetry have 
been also checked by inspection and confirmed for the transmission matrix (2.34) (see Ap- 
pendix B for details). In addition to the derivation of the transmission matrix eigenvalues 
via the BAE, there exist strong physical and algebraic arguments, as well as various checks 
regarding the validity of (2.34) as is manifest from the discussion above. It is also worth 
noting that this is the first time to our knowledge that such expressions are computed in 
the context of the XXX spin chain with a defect from the Bethe ansatz equations, and 
expressions (2.28), (2.34) are as far as we can tell novel. 



3 The anisotropic case: XXZ model 

Our investigation on the defects is carried on with the anisotropic XXZ spin chain. Based on 
the analysis of this model we shall be able to make contact with relevant results extracted in 
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the context of the sine-Gordon modeL We shall focus in the present study on the so called 
type-II defects [4]. Let us first introduce the general defect L- matrix [29]: 



^sinh (^^(A + iS' + £ 
sinh(z/i) 



sinh(z/i) 5* 
sinh ( yu(A — zS*^ + 



(3.1) 



where we define q = e*'^, measuring the anisotropy of the model as A = coshi/i. The bulk 
L-matrix is basically the i?-matrix of the XXZ model, that is the spin ^ representation of 
the latter expression, i.e. one simply implements in (3.1) the identifications (2.3). 

The L-matrix is the typical solution of the quadratic algebra associated to iXq{si2), which 
reads as: 



s+, s- 



^2S' — q ^•S"' 
q-q-^ 



(3.2) 



In the finite case, which will be considered here, the generators are represented by n x n 
(n = 25* + 1), matrices as: 

n n— 1 n— 1 

•S*^ = «fc Cfcfc , = Cfc Cfcfc+i , S = C*fc Cfc+ifc , (3.3) 



fc=l 



k=l 



k=l 



where we define 



ak = -{n + l-2k) 



Ck 



[kUn-k] 



Q 1 



q — q~ 
q-q- 



(3.4) 



Note that in [4] the L-matrix (3.1) is also used, but an infinite dimensional representation is 
employed; nevertheless, the structure of the defect matrices as well as the extracted physical 
quantities presented in [4] are similar to our expressions, as will be evident later in the 
text. The Hamiltonian of the model is given by (1.3) where L(0) = /i diag (^cosh ifi{S^ + 

i), coshi/i(-S^ + i)), 

hjj+i{0) = ^ Ujcr^^i + f^Jf^J+i + cosh(z/i)a|(Tj+i + cosh(z^)IjIj+i 



■'n+ln 



(0) 



sinh(z/i) 



^ + S^al^^) + cosh(y) sinh(z^5^X+i + sinh(y) cosh(z/i5;jI„+i. 



The generic defect matrix and the bulk L-matrix possess highest weight states, thus as in 
the isotropic case the typical algebraic Bethe variation may be applied and the corresponding 
Bethe ansatz equations are immediately obtained^ having the following standard form 



N 



ey{X, - 0) ef (A,) = - J] e2{Xi - A,) , y = 2S , 



(3.5) 



■^Thc BAE are valid for any S ^ real number. 
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with 6 being again the rapidity of the defect, and we also introduce the notation 

3mh(^(A + f)) 
- .inh(MA - f )) ■ 

We may now proceed with the derivation of the transmission matrix distinguishing two 
regimes, the repulsive and the attractive, depending on the value of the coupling constant. 
Bear in mind that we wish to compare our findings with similar results in the context of the 
sine-Gordon model. Therefore, it is useful to provide the relation between the sine-Gordon 
coupling constant and the anisotropy parameter fj, of the XXZ model (see also e.g. [30] 
for a more detailed discussion): 

= 8(7r — /i) , 47r < < Siri repulsive regime, 

= 8/i , < < 47r attractive regime. (3.7) 

Note that in the attractive regime the formulation of bounds states between solitons and 
anti-solitons of zero spin (topological charge), the so called "breathers" is allowed. 



3.1 The repulsive regime; the transmission matrix 

We shall first consider the repulsive regime, and derive the corresponding transmission ma- 
trix. In this regime the antiferromagnetic ground state is a filled Fermi sea with real strings 
(A^ odd), as in the isotropic case. For our purposes, it suffices to consider here the state with 
two holes. The density associated to this state is obtained in the thermodynamic limit by 
following the logic described in the previous section 

1 /"°° 1 ^ 

a(A) = ai(A) + -a,(A " ©) " y dX' a^iX - A') a(A') + ]^ Z ^^(^ - ^.O > (3-8) 

where the basic formula in the presence of two "holes" of rapidities Xj (2.14) was exploited. 
Note that a„(A) is defined again as in (2.13), whereas its Fourier transformation is given by 

^ ^ ^ sinh((z/ — n)^) n 

In the present subsection, we restrict ourselves to the values < y < 2v. Results on generic 
values of y are presented in Appendix C. 

As in the isotropic case, the density may be expressed in a compact form as 

a{X) = ao(A) + - ( ^ r,(A - A,.) + r,(A - 6)) . (3.10) 
fc=i 
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The Fourier transforms of the latter quantities have been exphcitly determined 

1 sinh((z/-2)f) sinh((z/-y)f) 



2cosh(f)' ""'^"^^ 2sinh((z/- l)l)cosh(f) ' 2 sinh((i. - l)f ) cosh(f ) 

(3.11) 

Relations (2.19) for the energy and momentum of the particle-like ("soliton") excitation are 
also valid here. 

In order to identify the scattering amplitude between two holes, we compare the expres- 
sion providing the density of the state with the quantization condition (2.21), with respect 
to the excitation of rapidity Ai. We are now in the position to compute the soliton-soliton 
scattering amplitude, as well as the transmission amplitude for the XXZ model using the 
expressions (2.22). Taking into account the following useful identity 

4 Jq X sinh x sinh fix n ^2~'~2~'~'^~'~2''' ^ ^ 



k=0 



we may easily reproduce the well known expression for the sine-Gordon soliton-soliton scat- 
tering (see also [31, 4]), 

fr r(z + 2(fc + l)7) r(z + 2A;7 + l) 
' Al r(z + (2A; + 1)7) r(z + (2A; + 1)7 + 1) 

^ r(-z + {2k + 1)7) n-z + {2k + 1)7 + 1) 

r(-z + 2(A; + l)7) r(-z + 2A;7 + l) ' ' 

where we now define 

Suitable configurations corresponding to soliton/anti-soliton states may be formulated, so 
that all the S'-matrix eigenvalues may be identified. We shall not give the details of such 
a derivation here, however we refer the interested reader to e.g. [30, 32] and references 
therein for a more detailed analysis. The S'-matrix, solution of the Yang-Baxter equation, is 
structurally similar to the "bare" i?-matrix, and is given by 



K\l) c(7) 
c(7) fc(A,7) 

V «(A,7)/ 



(3.15) 



where 

a(A, 7) = sin(7r7(zA + 1)) , /3(A, 7) = sin(z7r7A) , 0(7) = sin(7r7) . (3.16) 
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Note that the S'-matrix is essentially a renormalized i?-matrix, given that both the spectral 
parameter as well as the anisotropy parameter are renormalized, as is clear from the expres- 
sions above. The S'-matrix is in agreement with the resuhs on the sine-Gordon model (see 
also [31, 4, 30, 32]). 

Recalhng expressions (2.22) we can also extract the transmission amplitude. Based on 
the state described above, as well as the factorization argument, we find the first eigenvalue 
of the transmission matrix 

r(z + 7S + f + {2k + 1)7) r(2 - 75 - f + i2k + 1)7 + 1) 



T(A,7,^) = n 



T{z + ^S+^ + 2k-f) r(2 - 7S - I + 2{k + 1)7 + 1) 



' \i-z + 7^ + i + 2fc7) Tj-z - 7^ - 1 + 2(fc + 1)7 + 1) 
T{-z + ^S+^ + {2k+ 1)7) T{-z - 7^ - I + (2fc + 1)7 + 1) ' ^ ■ 

where z = iX'y and S = S — ^ is the "shifted" spin. As already mentioned, in order to derive 
the transmission amplitude it is sufficient to consider the state with one hole. In this case, 
the corresponding spin eigenvalue is computed explicitly and turns out to be 

S' = ^{S+h. (3.18) 
u — 1 2 

One immediately observes an overall renormalization factor ^j^^, which is equal to the ratio of 
the bare anisotropy parameter u over the renormalized one I' — l; the physical (renormalized) 
spin reduces then to the expected one: S^^^^ = S + ^. 

Alongside the S'-matrix (3.15), the transmission matrix T satisfies the quadratic algebra 
(2.33). We conclude that the transmission matrix may be cast then as 

jr^ ^ S)= Tjkl^S) f sin(7r7(iA + + \)) sin(7r7) S~ \ 

^ ' ^ sin(7r7(iA + ^ + i)) \ sin(7r7) 5+ sin(7r7(d - + i))^ ' ^ ' ' 

where , S"^ correspond to the spin S representation of ilg(sl2), where now q = e^'^^ . We 
should point out that the transmission matrix is essentially a renormalized defect matrix, as 
is manifest from the structure of T and L-matrices. Having said this, it is straightforward 
to verify that T is indeed a solution of the fundamental algebra (2.33). Moreover, based on 
the analysis presented in Appendix A, we have checked by inspection that the transmission 
matrix satisfies unitarity and crossing symmetry, hence its validity is completely confirmed. 
As in the isotropic case, appropriate string configurations suitably positioned with respect to 
Ai, provide the various eigenvalues of the transmission matrix. However, in the trigono- 
metric case this is a highly intricate task, and will be left for separate investigations. It 
is also worth noting that this is the first time as far as we know that such expressions are 
computed via the Bethe ansatz formulation in the trigonometric case. A detailed comparison 
with earlier results on the transmission matrix of the sine-Gordon model will be given in the 
next subsection. 
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3.2 The attractive regime 



3.2.1 The soliton transmission matrix 

We shall now focus on the attractive regime; in this regime bound states between solitons 
and anti-solitons, the so-called "breathers", exist. Thus the scattering between the breathers 
and the defect may be also investigated. As was shown in earlier studies, the ground state 
in the attractive regime consists of the so-called negative parity strings (see also e.g. [30] 
and references therein) 

A(-) = A + ^. (3.20) 

Zi LA] 

The BAE are modified then as follows 

M 

gy{K - 6) (A) = - n e2(A. - A,) , (3.21) 



where we define 



cosh(u(A + f )) 
9n{\) = .r^, .1 . (3.22) 

COSh(yU(A - yj) 



A generic state with two particle excitations (two holes in the filled Fermi sea of negative 
parity strings) is considered and the density associated to this state may be derived based 
on the standard formulation [20, 25]. It turns out that the derived state density is given by 
the following expression 



1 /"°° 1 ^ 

^(A) = &i(A) + -hy{\ - e) - / d\' a,{X - A') a(A') + $^a2(A - A,) , 



(3.23) 



where the formula (2.14) in the presence of two "holes" of rapidities Xj has been exploited. 
The Fourier transformation of a„ is given in (3.9), whereas the Fourier transform for bn is 
found to be 

sinh(^) 

M^) = --^y 0<n<., 
- , ^ sinh((n - 2z/)^) 

For the sake of simplicity, we consider here the case < y < z/. Results regarding generic 
values of y are presented in Appendix C. Similarly to the previous sections, the density is 
compactly expressed as 

a{\) = ao(A) + - ( 5^ r,(A - A,.) + r,(A - 6)) , (3.25) 
fc=i 
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whereas the Fourier transforms of the latter quantities are given by 

. . . 1 . . . sinh((z/-2)|) 

"""^"^^ 2cosh((z/-l)f) ' ""'^"^^ 2sinh(f)cosh((z/-l)f) ' 

sinh(?/|^) 

= 2sinh(f)cosh((z.-l)f) • ^^-^^^ 

Note also that relations (2.19) also hold for the energy and momentum of the particle- like 
excitation. 

Comparison of the expression providing the density of the state (3.25) with the quantiza- 
tion condition (2.21), with respect to the excitation of rapidity Ai, leads to the derivation of 
the soliton-soliton scattering as well as the transmission amplitude, given by (2.22). Taking 
into account the identity (3.12) we reproduce the well known expression for the sine-Gordon 
soliton-soliton scattering in the attractive regime (compare for instance with the notation 
used in [4]), which is given by (3.13), but now we define 

z = iX, 7 = z/-l. (3.27) 

Soliton anti-soliton configurations leading to all the S'-matrix eigenvalues may be identified 
[30, 32]. The S'-matrix, solution of the Yang-Baxter equation as well, is structurally similar 
to the "bare" i?-matrix, and is given as (3.15), where we now define 

a(A, 7) = sin 7r(iA + 7) , /3(A, 7) = siniTrA , 0(7) = sin 717. (3.28) 

In a similar fashion, we also identify the transmission matrix, which describes the inter- 
action between the soliton and the defect. Based on the state described above, as well as 
the factorization argument, we extract the first eigenvalue of the transfer matrix, which is 
expressed as 

. ^ fv r(z-e + 2(fc + i)7 + |)r(z + e + 2fc7 + |) 

^ J-jLr(5-e + (2A; + l)7 + |)r(i + e + (2fc + l)7+i) 

r(-^ - e + (2fc + 1)7 + I) T{-z + e + (2fc + 1)7 + 1) 

r(-i-^ + 2(/c + l)7 + i) r(-^ + e + 2A;7 + i) ' ^' ^ 

where 

z = tX, ^ = S+l- (3.30) 

The transmission matrix T satisfies the quadratic algebra (2.33) together with the S'-matrix 
derived above. We conclude that the transmission matrix T may be then cast as (set ^ = iu) 

TP/ ^(^^7,'^) fsm{7r-f{iu + + 1)) sin(7r7) \ 

sm(7r7(«M + S + ^)) \ sm(7r7) S^ sm(7r7(2M — S^ + 2)) / 
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The elements S^, form now the spin S = representation of ilg(sl2) with q™'^, which is 

an infinite dimensional representation. Notice that here we have used the fact that: 
s+- 

sin(7r7(m + | ^)) = ± sin(7r7(m + |)) or ± cos(7r7(m + |)) given that S* + | is an integer 

or half-integer. Computation of the spin of the state with one hole via the Bethe ansatz 
equations leads to S*^ = |, i.e. the renormalized spin is ^ (z/ is the renormalization factor 
in the attractive regime), which is basically the spin of the hole, denoting that the defect 
spin is effectively zero, confirming the relevant result {S = 0) through the derivation of the 
transmission matrix. 

To efficiently compare with earlier relevant results from the Sine-Gordon model [4], it is 
convenient to introduce some useful notation; first we shift A such that 

iX-^iX + A, (3.32) 

where A is an arbitrary constant. Also define the following quantities: 

Zj = -tX-^7]j, r/i = -(A + 0, r/2 = -(A-0, J e {1, 2}. 

TT 7 7 

(3.33) 



The transmission factor is then expressed in terms of Zi , Z2 as 



TT 

where the overall physical factor pd is given by 



sin7r(z2 + |) 

T{zi, Z2) = pd{zi, Z2) , (3.34) 



, , p.l .f^r(^i + (2A;-l)7 + i)r(^2 + (2A:-l)7+i 



'2 '2 '^11 T{z, + 2k^ + l)T{z2 + 2kj + l) 

r(-^i + 2fc7 + i) r{-z2 + 2k^ + l) 

Y(-z, + (2k - 1)7 + i) T{-Z2 + (2A; - 1)7 + I) ■ ^ ' ^ 

The latter expression is identical to the one extracted in [4] for the type-II defects, taken 
into account the identifications (3.33). Having determined the soliton transmission matrix 
we now proceed with the derivation of the breather transmission matrix. 



3.2.2 The breather transmission amplitude 

The breathers are in general identified within the Bethe ansatz frame by suitable string 
configurations. To fully describe this scattering process for the breathers it is necessary to 
take into consideration two sets of Bethe ansatz equations; the first set describes the negative 
parity one-strings, while the second one describes the breather itself. The second set of BAE 
is necessary in order to derive the energy and momentum of the breather, and also compare 
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with the quantization condition with respect to the breather (for more details on breathers 
and their interactions we refer the interested breather to [30] and references therein). 

A state with two hght breathers with rapidities Ai, A2 will be considered. We shall basi- 
cally deal with the lightest breathers for simplicity; a generalization of the results concerning 
higher breathers is then straightforward [30], and is given at the end of the subsection. The 
lightest breather is described by one positive parity (real) string with rapidity A,, then the 
BAE for the state with two breathers are expressed as : 

M 2 

gyiX, - 6) 5f (A,) = -ne2(A. - A,) l[g2iX^ - A,). (3.36) 
There is a second set of BAE describing the breather with rapidity Ai 

M 

e,(Ai - 6) ef (Ai) = - n^?2(Ai - A,) e^CXi - X2). (3.37) 

i=i 

As already mentioned the second set is necessary for our purposed here, given that it facili- 
tates the computation of the energy and momentum of the breather as well as the formulation 
of the corresponding quantization condition. 

From the first set of BAE (3.36) the following density regarding the negative parity 
strings arises, 

1 ^ 

(A) = MX) + - (b(A - 0) + ^ R(A - A,)) , (3.38) 



2 

a 

/V V ^ 

i=i 

where we define the Fourier transforms of R, B as 



cosh(f ) . sinh(f ) 

^^""^ ~ "cosh((i.-l)f)' ^^^^ - 2sinh(f)cosh((z.-l)|) ^^■^^> 

The second set (3.37) leads to the density describing the breather 

2 

(t(A) = ao(A) + {t,{X - 9) + Y,n{X - A,)), (3.40) 

i=i 

where the corresponding Fourier transforms read as 

. ^ ^_ cosh((^-2)|) cosh((^-3)f) cosh(0^-y-l)f) 

- cosh((^.-l)|)' ""^^^ - "cosh((^-l)f)' ^"^^^ - cosh((z.-l) 

Moreover, it may be shown as in the soliton case that 



' (3.41) 



ao(A)=e(A), ^lA) = i-^ (3.42) 
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where e and p are the energy and momentum of the hghtest breather respectively. 
Similarly, a quantization condition for the breather may be formulated 

^ivix^)Ng^^^^ A2, e) - 1) |Ai, A2) = 0, (3.43) 

due to the factorization of the scattering process S'(Ai,A2,0) = 5'^^'^''(Ai, A2) ^b^^(^i7 0)- 
Comparison of the latter formula with a leads to the expressions for the breather scattering 
amplitude as well as the corresponding breather transmission amplitude: 

5p)(A)=exp[- r —e-'^^h{^)\. t/'^(A) = exp [ - H —e-'^\{uj)\ (3.44) 

L J~oo ^ -I L J-00 ^ -I 

A = Ai - A2, A = Ai - e. 

Bearing in mind (3.41), (3.44) as well as the useful identity (2.23) we conclude that the 
breather scattering amplitude is given by the following hyperbolic ratios, 

_ smh(l-g) .inh(| H-fH-f) 
smh(| + i|) smh(f-i|-f)' 

where we define: 9 = The scattering amplitude coincides of course with the familiar 
sine-Gordon breather quantity, and for the lightest breather, which we consider here this 
corresponds to the scattering amplitude of the scalar sinh-Gordon field [31, 4, 30]. 

Similarly, through (3.44) the breather transmission amplitude may be derived as 



sinh(^ + ^) sinh( 



zUl _|_ HE^ 
4 ^ —"V 2 ~'~ 4 



9 = the constants rji have been defined in (3.33), and it is clear that our expression (3.46) 
for the breather transmission amplitude coincides with the one identified in [4]. 

The results on the scattering and transmission amplitudes may be generalized for higher 
n-breathers, which are represented by n-positive parity strings with real centers Xj. More 
precisely, it is straightforward to see (see e.g. [30]) that the scattering between two generic 
rii, n2 breathers may be expressed as: 



ni n2 



St''''\>^) = \[ n ^i'''^(A + ^(ni-n2-2/, + 2/2)), (3.47) 



il=l l2 = l 

while the transmission amplitude of an n-breather is 



n 

Tr\X) = X{Tl^'^{\ + \{n + l-2l)). (3.48) 

1=1 
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The spin of the one n-breather state in the presence of the defect can be also evaluated, 
and turns out to be = 0. The spin of the breather is as expected zero, recall that in the 
attractive regime as pointed out in the previous subsection the "effective" defect spin is zero. 
With this we conclude our analysis on the transmission matrices for the kinks (solitons) and 
breathers in the (an) isotropic Heisenberg models. 

4 Discussion 

We have studied in the present article one-dimensional (an) isotropic Heisenberg chains in 
the presence of a single impurity. We have recalled the construction of the kink S'-matrix 
and have produced the transmission amplitudes through the Bethe ansatz equations. We 
have also been able to derive the breather's transmission matrix in the attractive regime of 
the XXZ chain. Our findings in the attractive regime of the XXZ model coincide with earlier 
results obtained in the context of the sine-Gordon model [4], implying that the picture is 
consistent. Our findings in the XXX C3iSG clS well as in the repulsive regime of the XXZ 
model are novel as far as we can tell and involve finite representations of the and iig(sl2). 
As a further check, comparison with relevant results at the classical level should be made, 
especially with studies in the Landau-Lifshitz [19], and sine-Gordon models [17]. 

Some comments on future directions are in order here. We have based our formulation 
here on finite representations of 5(2 , ilg(sl2). More precisely, we have restricted our analysis 
on defect matrices, that correspond to finite representations, thus the algebraic Bethe ansatz 
formulation may be applied. The chosen representations possess high weight states, and 
therefore the algebraic Bethe ansatz variation can be successfully applied. However, infinite 
dimensional representations of 5(2, Ug(s[2) for the defect matrix, can be also considered. In 
this case local gauge (Darboux) transformations in the spirit described in [33] should be 
employed in order to extract the associated spectrum and BAE. 

Moreover, there exist several open issues to be resolved; it is an intriguing task to de- 
termine the string configurations as well as the positions of the real centers of the n-strings 
within the XXZ model, in order to extract all the eigenvalues of the trigonometric transmis- 
sion matrix from the BAE. A natural generalization would also be to extend our analysis in 
the case of higher rank (deformed) algebras. Finally, depending on the values of the coupling 
constant it is possible to consider the formation of bound states between the particle-like 
excitations and the defect. This analysis may be achieved via the investigation of the poles 
appearing in the overall physical factor of the transmission matrix. All the above are signif- 
icant issues, which hopefully will be addressed in the near future. 
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A Eigenvalues & eigenstates of the defect matrix 



In the present appendix we identify the eigenvalues and the corresponding eigenstates of the 
defect matrix. Note that the spin operator shares essentially the same set of eigenstates. We 
begin with writing the defect matrix as 

n-l 



m=l m=l 

L{X) G End^C^ ® ^^[A] n = 2S + 1. We consider the following ansatz for the eigenstates: 

iV^fc) = /fc e? ® 4=^1 + /fc+i ® 4"^ = 1, ■ ■ ■ , - 1 , (A.l) 

where e^'^'' is a rf- dimensional column vector with zero elements everywhere, and the unit at 
the k-th position. 

Acting with L on the vector, according to the rule e^f^j ef^ = 6hc (fa \ we find 

L(A) IV^o) = (A + l^o), ^(A) l^n) = (A + y) l^n) , (A.2) 

and 

L(A) l^fc) = . . . = /fc4'^®4+i(A+i+^fc?/fc) + /fc+i4'^®4"^(^^:+Cfcyfc') = en,fc l^fc) , (A.3) 
where we have defined yk = A^. In order for lipk) to be an eigenvector we require the 

Jk 

following relation to hold 

Ak+i + BkVk = "Pfe + CkVk^ ' 

or equivalently, 

Bk yl + yk {Ak+i - Vk) - = . 

Solving this algebraic equation provides us directly with the eigenvalues of the defect oper- 
ator, which would have the generic expression 

Rational case: Let us first define the following functions in the rational case 

Ak = \ + iak + -, Bk=Ck=iCk, Vk = \ - iak + - . (A.4) 

Substituting the corresponding functions for the rational case leads to the following values 
of yk 



Vk 



(1) I k 1 



n-^^ ,(2) 
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implying that the corresponding eigenvalues are 

(1,2) _ , 

It is interesting that the eigenvalues are fc- independent, thus the generic spin-S* defect matrix 
possesses only two distinct eigenvalues, each one being n-fold degenerate. 



Trigonometric case: The relevant functions in this case are: 



Bk=Ck = ^(g - q ^)Ck, 



2 



(A.5) 



After some calculations, we find that the eigenvalues in the trigonometric case are much 
more complicated, and given by the following expression 



,(1,2) 



1 



cos[|(n— 2/c)] sinh(A/i)±- —l—cos[jj,{2k—n)]+2 cos(n/i)+cosh2A/i(— l+cos[/i(n— 2/;;)]) 

^ (A.6) ■ 



Spin eigenvalues: We may also compute the total spin associated to the defect matrix 

= y ® In + I2 ® S', (A.7) 

where 5*^ is given by the n x n matrix (2.4) corresponding to the spin 5* representation, and 
Irn are the m-dimensional unit matrices. As noted already, the states found above (A.l) are 
also the eigenstates. The spin eigenvalue problem reads then as 

|z^o)= (5+^)1^0), 

m ={s + ^-k) liJk), A; G {1, . . . , n - 1} 

l^n) = ( - 5 - ^) \^n) . (A.8) 

We end up with 2n eigenstates and eigenvalues. Each eigenvalue S — k + ^, A;G|l,...,n— 1| 
is 2-fold degenerate, so that in total there are n + 1 distinct eigenvalues, as expected from 
the spin summation rules, which take the familiar values: E < — S* — |,...,5' + |k 



B Unitarity, crossing symmetry & Casimir operators 

In this appendix we present the unitarity and crossing symmetry properties of the transmis- 
sion matrix, and confirm the findings of the present investigation by exploiting these basic 
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requirements together with the fact that the transmission matrix satisfies the quadratic al- 
gebra (2.33). The unitarity and crossing symmetry are given by the following expressions 
respectively 

Ti2(A) Ti2(-A) =1, 

T'lU-A + i) Ml T%{X + i) Ml = I , (B.l) 

where ti denotes transposition on the first vector space. In both the isotropic and anisotropic 
case in the principal gradation, which are considered here, M = I. 



Rational case: Let us first consider the XXX isotropic case. Recall the expression found 
for the transmission matrix associated to the spin S representation: 

T(A,^) = ^^^^^^ , M(A, S) 

where T(A) is defined in (2.28). 

The requirement for unitarity leads to: 

T(A,^) T(-A,^) = 1, (B.3) 

whereas the crossing symmetry condition provides: 

X . X (^A + 5 + i) (-zA + 5 + i) , , 

T A + S) T -A + 2, S) ) = I. B.4 

[iX + b - i) [-iX + 5 - ^) 

To obtain the relations above we have used the following 

Mi2(A, S) Mi2(-A, S) = M{\{X + S) M\\{-X + z,S) = 

(A2 + C)I= (zA + 5 + i)(-zA + 5 + i)I (B.5) 

C is the Casimir operator of SI2, and one may easily show for the spin S representation that 

C = (S^^ + Us-, sU + ] = ^^^±ll, (B.6) 



2 I ' J 4 4 

The latter is valid even for any generic real S. A discussion on the spin S (any real) 
representation of expressed in terms of differential operators may be found in e.g. [23] 
and references therein. It is easily confirmed by inspection that T defined in (2.28) satisfies 
relations (B.3), (B.4) emanating from the basic properties. 
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Trigonometric case: Similarly for the anisostropic XXZ chain recall the transmission 
matrix associated to the spin S representation of the deformed algebra 

T(A,/i,^) = -T^nr?TTTT ^) 
sm(/i(«A + 6 + 2)) 

' y sin(/i) 5+ sm{i2{tX- + l))J ^ ' 

H = 7r7, T(A,5) defined in (3.17), (3.31). 

As in the isotropic case the unitarity and crossing-unitarity conditions lead to: 
T(A,/i,^) T(-A,/i,5) = I, 

, sin(u(2A + S" + i)) sm(u(—i\ + S+^)) 
T(A + z,/i,g)T(-A + z,/i,g) . r^.,^^ . ..^^ =1 (B.8) 

sm(/i(«A + — 2)) sm(/i(— zA + — 2)) 

Again confirmation of the findings regarding the anisotropic case is easily done by inspection. 
The latter relation were obtained using the fact that: 

Mi2(A, fi, S) Mi2(-A, /i, S) = M*^2(A + ^, /i, S) M%{-\ + z, /i, S) = 
- ^ cos{2fiiX) + = sin(/i(zA + S + ^)) sin(/i(-zA + 5 + ^))I (B.9) 

Cg is the associated q-deformed Casimir operator of the \iq{s{2) algebra (g = e*^), and one 
may easily show for the spin S representation that 

= gg25^+g-lg-2^?^ + (g-g-l)2 5" 5+ = gg-2S^ +^"1^25^ + (g_g-l)2 5+ 5- = 2 C0s(/i(2S+l)). 

(B.IO) 

The results for both rational and trigonometric cases are valid for any real values S of the 
representation. In this case the representation may be expressed in terms of differential 
operators (see e.g. [23] and references therein). 



C The XXZ model: generic values of the spin 

In the main text we have computed the transmission matrices given that the spin parameter 
takes values in a restricted interval. We generalize here our results for generic spin values. 
It is clear that these generalizations are due to the periodicity of the involved trigonometric 
functions. 



Repulsive regime: Let us introduce the Fourier transform of ay for generic values of 
y = 2S: 

sinh ('((2m+ l)z/-?/)f) 
ay{uj) = ^ . , 2mp < y < 2{m+l)p. (C.l) 

Slllil ^Z/ 2 ) 
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The effect of the latter generahzation has the following effects in our computations in the 
repulsive regime of the XXZ model: the new "shifted" spin becomes 

S = S-m-]^ (C.2) 

whereas the expression for the transmission matrix is again given by (3.19), where 

" T{z + ^S -m+l + {2k + 1)7) r(5 - 75 + m - I + (2/c + 1)7 + 1) 



T(A,7,^)= = n 



T{z + -iS -m + l + 2k^) T{z--fS + m-^ + 2{k+ 1)7 + 1) 
Tj-z + 7g - m + I + 2A:7) Tj-z - 7^ + m - | + 2(/c + 1)7 + 1) 
^ T{-z + -m + ^ + {2k + 1)7) T{-z - 7^ + m - | + (2A; + 1)7 + 1) 

z = iX'-f. 



y 



Attractive regime: First we introduce the generalized Fourier transform for h 

« , , sinhffw — 2mz/)^) , ^ , ^ 

byiu) = ^l^^^^^^ ^^ < y < {m+l)u. (C.3) 

Then the transmission matrix in the attractive regime is given by (3.31), where 

T{z-^ + m(7 + 1) + 2{k + 1)7 + ^) r(z + e - m(7 + 1) + 2A;7 + i) 



T(A,7,^) = n 



X 



,=0 -^ + m{^ + l) + {2k + 1)7 + 1) r(£ + ^ - m(7 + 1) + {2k + 1)7 + \) 
V{-z - e + m(7 + 1) + {2k + 1)7 + \) V{-z + e - m(7 + 1) + {2k + 1)7 + \) 
T{-z - e + m(7 + 1) + 2{k + 1)7 + \) T{-z + ^ - ^(7 + 1) + 2A;7 + \) 



2' 

where 



z = i\, ^ = 5 + |, S = m. (C.4) 
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